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Abstract. Given a graph G, the number of nowhere-zero Zq-flows 4>g(o) is known to be a polynomial in q. We 
extend the definition of nowhere-zero Z q -flows to simplicial complexes A of dimension greater than one, and prove 
the polynomiality of the corresponding function <f>A(q) for certain q and certain subclasses of simplicial complexes. 

Resume. Efant donne une graphe G, on est connu que le nombre de Z 9 -flots non-nuls G (q) est un polyn~ome dans 
q. Nous etendons la definition de Z g -flots non-nuls pour incluir des complexes simpliciaux de dimension plus grande 
qu'un, et on montre que le nombre est aussi un polyn'ome de la fonction correspondante pour certain valeurs de q et 
de certaines sous-classes de complexes simpliciaux. 
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1 Introduction 

Flows on graphs were first studied by Tutte ( [Tut47, Tut48]; see also [BSST75] and [Tut76]) in the 
context of Kirchoff's electrical circuit laws. Since this time, a great deal of work (see, for instance, 
[Jae88], [Sey95], [Whi92, Ch. 6, Sec. 3.B, 3.C], or [BZ06]) has been done, and applications found in 
network and information theory, optimization, and other fields. 

To define a flow on a graph G, we first give an initial orientation to its edges (this orientation is arbitrary, 
but fixed). Then, a Z q -flow on G is an assignment of values from Z 9 to each edge such that modulo q, the 
sum of values entering each node is equal to the sum of values leaving it. If none of the edges receives 
zero weight, the Z g -flow is nowhere-zero. In terms of the signed incidence matrix M of G, a Z g -flow is 
an element of the kernel of M mod q. 

In this paper, we extend the idea of flows on graphs to flows on simplicial complexes (namely, elements 
of the kernel of the boundary map mod q — we give detailed definitions below), and explore the work 
done on graph flows in the context of simplicial complexes. We consider this paper as the starting point 
of a systematic study of flows on simplicial complexes. Extending Tutte's polynomiality result for the 
enumeration of nowhere-zero flows on graphs, our first main result is the following: 

Theorem 1 Let q be a sufficiently large prime number, and let A be a simplicial complex of dimension 
d. Then the number 4>a(q) of nowhere-zero Z q -flows on A is a polynomial in q of degree /3d(A) = 
dimQ(iJ d (A, 
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In other words, there exists a polynomial p(x) that coincides with 4>&(q) at sufficiently large primes. 
In fact, 4>\(q) coincides with p(x) at many more evaluations: 

Theorem 2 The number 4>a(q) of nowhere-zero "Lq-flows on A is a quasipolynomial in q. Furthermore, 
there exists a polynomial p(x) such that <j)/±(k) = p(k) for all integers k that are relatively prime to the 
period of 4>\{q). 

(Again, we give detailed definitions, e.g., of quasipolynomials and their periods, below.) For the case 
where the simplicial complex triangulates a manifold, we prove more. 

Proposition 3 Let Abe a triangulation of a manifold. Then 

if IS. has boundary, 

q — 1 if A is without boundary, Z-orientable, 

if A is without boundary, non-Z-orientable, and q ^ 2, 

1 if A is without boundary, non-lL-orientable, and q = 2. 

In Section 2, we provide background on simplicial complexes and boundary operators, and define 
flows on simplicial complexes. We also discuss certain details of the boundary matrix of pure simplicial 
complexes. In Section 3, we prove Theorems 1 and 2. In Section 4, we examine the specific case of 
triangulated manifolds and prove Proposition 3. In Section 5, we present several open problems. 



<M(<?) = < 



2 Flows on Simplicial Complexes 

Recall that an (abstract) simplicial complex A on a vertex set V is a set of subsets of V. These subsets 
are called the faces of A, and we require that 

(1) for all v 6 V, {v} 6 A, and 

(2) for all F G A, if G C F, then G e A. 

The dimension of a face F is dim(F) = |F| — 1, and the dimension of A is dim(A) = max{dim(_F) : 
F £ A}. A simplicial complex is pure if all maximal faces (that is, faces that are not properly contained in 
any other face) have the same cardinality. In this case, a maximal face is called a facet, and a ridge is a face 
of codimension one. For more background on simplicial complexes, see [Sta96, Ch. 0.3] or [Hat02, Ch. 
2]. 

Let A be a pure simplicial complex of dimension d with vertex set V = {v , ...,«„}. Assign an 
ordering to V so that v < vi < ■ ■ ■ < v n . Any r-dimensional face of A can be written (in terms of this 
ordering) as [vi ■ ■ ■ Vi r }. Let d be the boundary map on the simplicial chains of A, defined by 

r 

9[v io ■ ■ ■ v lr ] = K ■ • • ■ ■ ■ v ir ]. 

3=0 

Define the boundary matrix dA of A to be the matrix whose rows correspond to the ridges and columns 
to the facets of A, and whose entries are ±1 or 0, depending on the sign of the ridge in the boundary of the 
facet. In the case of a graph, the facets are the edges, and the ridges are the vertices, and, under the natural 
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ordering of the vertices, the boundary matrix is identical to the signed incidence matrix of the graph. For 
more background on the boundary map and boundary matrices, see [Hat02, Ch. 2]. 

Now, since a flow on a graph (one-dimensional simplicial complex) is an element of the kernel mod 
q of the signed incidence (boundary) matrix, we have a natural way to extend the notion of a flow on a 
graph to a flow on a simplicial complex. 

Definition A Jj q -flow on a pure simplicial complex A is an element of the kernel of d A mod q. A 
nowhere-zero Z q -flow is a Z g -flow with no entries equal to zero modq. 

As an example, consider the surface of a tetrahedron, with vertices V = {1, 2, 3, 4}, with the natural 
ordering. Its boundary matrix d A is: 





124 


134 


234 


123 


14 


-1 


-1 








24 


1 





-1 





34 





1 


1 





12 


1 








1 


13 





1 





-1 


23 








1 


1 



One instance of a nowhere-zero Z^-flow on A is (3:123, X124, X134, X234) 71 = (1, q — 1, 1, q — 1) T . 

For a second example, consider a pure simplicial complex A with a vertex v that is contained in every 
facet, i.e., A is a cone over v. 

Lemma 4 Let A be a pure simplicial complex that is a cone over a vertex v. Then A has no nontrivial 

Z q -fl0WS. 

Proof: Order the vertices v , v\, . . . , v n so that v = v n is the largest. Recall that the link of v in A is 
defined as 

lk A (w) := {F e A : v (£ F and {v} U F G A} , 

and the deletion ofv in A is 

A - v := {F e A : v g F} . 
Now consider the boundary matrix dA of A, with the rows and columns ordered in the following way: 

(a) Let all ridges containing the vertex v n come first, and let them be ordered lexicographically. 

(b) Let all ridges in lkA(vn) be next, ordered lexicographically. 

(c) Order the remaining ridges lexicographically. 

(d) Let all facets containing v n come first in the columns, and order them lexicographically. 

(e) Order the remaining facets (those that do not contain v n ) lexicographically. 

It is not hard to see that under these conditions the boundary matrix takes the following form: 
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9(lk A K)) 





±1 














Since A is a cone over v n , the two regions on the right-hand side of the matrix are nonexistent, and ±1 
extends through all columns. Therefore, we have a series of rows with precisely one nonzero entry. □ 



3 The Number of Nowhere-Zero Z^-Flows on a Simplicial Com- 
plex 

3. 1 Tutte polynomials of matroids and proof of Theorem 1 

To study the question of counting nowhere-zero flows, we recall the notion of a matroid and refer to, 
e.g., [Oxl92], for relevant background. Associated with every matroid M = (E,X), where E is the 
ground set and I the collection of independent sets of M, is the Tutte polynomial 

T M (x, *)=$>- iY iM) - r{s) (y 

SCE 

where r(S) denotes the rank of S. The Tutte polynomial satisfies the following recursive relation: 

{T M / e (x, y) + T M _ e (x, y) if e is neither a loop nor a coloop, 
xT M/e (x,y) if e is a coloop, 

yT M/e (x,y) if e is a loop. 

The Tutte polynomial has been generalized [OW79] by allowing the recursive definition to include 
coefficients in the case when e is neither a loop nor a coloop. We say that a function on the class of 
matroids is a generalized Tutte-Grothendieck invariant if it is a function / from the class of matroids to a 
field F such that for all matroids M and N, f(M) = f(N) whenever M = N, and for all e e E{M), 

J af(M — e) + rf(M/e) if is e neither a loop nor a coloop, 
~ \f(M(e))f(M - e) otherwise, 

where M(e) is the matroid consisting of the single element e, and a and r are nonzero elements of F. Let 
I be the matroid consisting of a single coloop and L be the matroid consisting of a single loop. We will 
use the following fact [OW79, Theorem 6.2.6]: 

Theorem 5 Let a and r be nonzero elements of a field F. Then there is a unique function t' from the class 
of matroids into the polynomial ring F[x, y] having the following properties: 
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(i) t'j(x,y)=x and t' L {x,y)=y. 

(ii) If e is an element of the matroid M and e is neither a loop nor a coloop, then 

t r M{x,y) = at' M _ e {x,y) +rt' M/e {x,y). 

(iii) If e is a loop or a coloop of the matroid M, then t' M (x, y) = t' M ,, (x, y) t' M _ e (x, y). 

Furthermore, this function t' is given by t' M {x, y) = a^ E ^ r ^ E ^T rt ^ E ^ Tm(^, ^) • 

We refer the reader to [Tut84], [Wel76], or [Whi92, Ch. 6, Sec. 1,2], for more background on and 
applications of the (generalized) Tutte polynomial. 

Finally, recall that any matrix may be realized as a matroid by taking the ground set to be the list of 
columns, and the independent sets to be the linearly independent subsets of columns. If E is the ground 

I £ I 

set of this matroid, and y G Z q , then we define the support of y as supp(y) := {e G E : y e ^ 0}. 

Proof of Theorem 1: Let A be a pure simplicial complex, and d A be the boundary matrix associated 
with A. Let M be the matroid given by the columns of d A, denoted as the ground set E. 
First, we claim that 

T M (0, 1 - q) = \{y G kcr(M) mod q : supp(y) = E}\ 
for any prime q that is sufficiently large. We show this by proving that the function 

9M(q) ■= \{y G kcr(M) mod q : supp(y) = E}\ 

is a generalized Tutte-Grothendieck invariant with a = — 1 and r = 1. The matrix for the single coloop 
/ is a single (linearly independent) vector, thus gi(q) = 0. The matrix for the single loop L is the zero 
vector, so = q — L Thus 5m(?) is well defined if \E\ = 1. 

Assume g is well-defined for \E\ < n, and let \E\ = n. Let e G E, and suppose that e is neither a 
loop nor a coloop. In the case of a matroid corresponding to a vector configuration (as is the case with 
the boundary matrix), contraction models quotients: let V be the vector space given by E (the columns of 
the boundary matrix), and let n : V — > V/V e be the canonical quotient map, where V e is the vector space 
with basis e. Then, the contracted matroid M/e is the matroid associated with the vector configuration 
{n{v)} ve E-{ e } m the quotient space V/V e . Let 

W := {y e kcr(M) : supp(y) = E}, 

X := {y G ker(M - e) : supp(y) = E - {e}}, and 

Z := {y G ker(M/e) : supp(y) = E — {e}}. 

We see that W n X = 0. By linearity of it, we have that W U X C Z. Now, suppose we have y G Z. 
This is a linear combination of all vectors in E — {e} which equals a scalar a G 7L q times e (that is, a 
linear combination equivalent to zero in the quotient space). If a ^ 0, then (y, a) G W. If a = 0, y G X. 
Therefore 

9ilf(?) = 9M/e(q) - 9M-e(q) ■ 
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If e is a loop (a column of zeros), then e may be assigned any of the q — 1 nonzero values of Z q . Thus, 
for loops 

9M(q) = (q- l)g M -e(q) 
= 9L(q)gM-e(q) ■ 

Finally, suppose e is a coloop. Then every maximally independent set of columns contains e, and we 
may, using row operations, rewrite M so that the column corresponding to e has precisely one nonzero 
element, and the row containing this element also has precisely one nonzero element. Therefore y e = 
for all y <E ker(M), and 

9M(q) = • gM-e(q) 
= 9i{q) 9M-e(q) ■ 

Thus gu (q) is well-defined and a generalized Tutte-Grothendieck invariant. We showed the cases for 
= 1 above, so we see that by Theorem 5 

9M(q) = t' M (0, q-l) = (-1)I £ |-^T M (0, 1 - q) 

and so 

T M (0, l-q) = \{y€ ker(M) mod q : supp(y) = E}\ . 

It follows that the number of nowhere-zero Z^-flows on a simplicial complex A is equal to Tsa {0,1— q) 
and hence is a polynomial in q. Using the definition of the Tutte polynomial, we see that the degree of 
this polynomial, in terms of the matroid, is \E\ — r{M). From linear algebra, we know that 

\E\ = dim(rank(M)) + dim(ker(M)) = r(M) + dim(kcr(M)) , 

so \E\ — r(M) = dim(ker(M)). But dim(ker(M)) is just the dimension of the top rational homology of 
A, dim Q (if d (A; Q)). By definition, this is /3 d (A), where d = dim(A). □ 

Remark. We require that q be prime as otherwise V and V e would be modules, rather than vector spaces. 
Requiring that q be sufficiently large ensures that the matrices reduce correctly over F^; that is, we require 
that the linear (in)dependencies are the same over Q as they are over ¥ q . For a simplicial complex of 
dimension d, a sufficient bound would be {d + 1)~^~, though this is not necessarily tight. 

3.2 Ehrhart quasipolynomials and proof of Theorem 2 

Let V C R d be a rational polytope, i.e., the convex hull of finitely many points in Q d . Then Ehrhart's 
theorem [Ehr62, BR07] says that the lattice-point counting function 

ehrp(t) :=#(fPnZ d ) 

is a quasipolynomial in the integer variable t, i.e., there exist polynomials p n {t),pi{t), . . . ,p/--i{t) such 
that 

ehi-p (i) = pj (t) if t = j mod k . 
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The minimal such k is the period of ehr-p(f) and the polynomials Po{t),pi{t), . . . ,pk-i(t) are its con- 
stituents. Since rational poly topes are precisely sets of the form {x <E R d : Ax < 6} for some integral 
matrix A and integral vector b, it is a short step to realize that the flow counting function </>a(<?) is a 
quasipolynomial in q: it counts the integer lattice points x = (x\, X2, ■ ■ ■ , Xd) (where d is the number of 
facets of A) that satisfy 

< Xj < q and dA(x) = mq for some m 6 Z . 

(Note that we only need to consider finitely many m.) So a (9) is a sum of Ehrhart quasipolynomials and 
thus a quasipolynomial in q We remark that a similar setup was used in [BS 1 lb] to study flow polynomials 
of graphs from an Ehrhartian perspective. 

Proof of Theorem 2: We just showed that a (<?) is a quasipolynomial, say of period p. By Dirichlet's 
theorem, there exist infinitely many primes of the form j + kp for gcd(j,j?) = 1 and k £ 1>q, and thus 
4>A (j + kp) agrees with the polynomial found in Theorem 1 for those j with gcd(j, p) = 1. □ 

3. 3 Quasipolynomiality of the Flow Function 

Given that the number of nowhere-zero Z q -flows on a graph is a polynomial in q, it is natural to ask 
whether the same is true for all simplicial complexes; however, this is not always the case. Consider the 
case of a Klein bottle, K. We have the following top homologies for K: 



H 2 (K,Z q ) = 



if q is odd, and 
Z 2 if q is even. 



Therefore, the number of nowhere-zero Z g -flows on K is given by a quasipolynomial of period 2: 



if q is odd, and 

1 if q is even. 



3.4 Calculating a Flow Polynomial: An Example 

Consider the triangular bipyramid A with vertex set {0, 1, 2, 3, 4} and facets {012,013, 023, 123, 124, 134, 
234}. Then dA is 



The kernel is generated by 





012 


013 


023 


123 


124 


134 


234 


01 


1 


1 

















02 


-1 





1 














03 





-1 


-1 














12 


1 








1 


1 








13 





1 





-1 





1 





23 








1 


1 








1 


14 














-1 


-1 





24 














1 





-1 


34 

















1 


1 



{(1, -1, 1, -1, 0, 0, 0) T , (0, 0, 0, 1, -1, 1, -if}. 
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Thus, any nowhere-zero Z^-flow will have the form (a, —a, a, —a + 6, —6, &, — fe) T where a,b G Zfe 
and a b. We have fc — 1 non-zero choices for a, and fc — 2 non-zero choices for 6. We see that 
4>/\{k) = (fc — l)(fc — 2). For primes fc > 1 (1 being the maximum over all absolute values of all 
subdeterminants of d A), we may also compute 4>A(k) = Tqa{0, 1 — fc) as 

SCE 

= fc 2 - 3fc + 2 
= (fc-l)(fc-2). 

4 Flows on Triangulations of Manifolds 

Proof of Proposition 3: Consider a pure simplicial complex that is the triangulation of a connected 
manifold. 

First, if a manifold has boundary, then there exists at least one ridge that belongs to only one facet. 
Since this corresponds to a row in its boundary matrix with precisely one nonzero entry, any vector in 
the kernel must have a zero in the coordinate corresponding to the facet containing this ridge. Therefore, 
manifolds with boundary do not admit nowhere-zero flows. 

If a triangulated manifold M is without boundary, then every ridge belongs to precisely two facets. 
This corresponds to every row of dM having exactly two nonzero entries. Therefore, since our manifold 
is connected, in a valid flow the assigned value of any facet is equal or opposite mod q to the assigned 
value of any other facet. It follows that every flow that is somewhere zero must in fact be trivial. 

It is known (see, for instance [Hat02, Chapter 3.3, Theorem 3.26]) that the top homology of a closed, 
connected, and Z-orientable manifold of dimension n is 

H n (M,Z) =Z. 

In terms of boundary matrices, this means that the kernel of the boundary matrix has rank one, since there 
are no simplices of dimension higher than n. By our comment above, we see that this element (a linear 
combination of the facets), must be nowhere-zero, and all entries are ±a G Z. It is then easy to see that 
triangulations of closed, connected, orientable manifolds have precisely q — 1 nowhere-zero Z g flows. 

If M is non-orientable, connected, and of dimension n, then we know (again, see [Hat02]) that the top 
homology is 

H n (M,Z) = 0, 

so there are no nontrivial flows on triangulations of such manifolds, except when q = 2. As every manifold 
is orientable over Z2, and since every row of the boundary matrix has precisely two nonzero entries, the 
all-ones vector is a nowhere-zero Z 2 -flow on M (in fact, the only one). □ 

5 Open Questions 

We have shown above that not all simplicial complexes have flow functions that are polynomials. How- 
ever, it is still natural to ask for which complexes 4>a(q) is a polynomial. Early experiments indicate 
this should include simplicial complexes that are convex-ear decomposable (see [Cha97]), which include 



Flows on Simplicial Complexes 



9 



complexes such as matroid, coloring, and a large subclass of unipolar complexes. We intend to use the 
structure of these complexes to further specify their specific flow functions. 

Many graph polynomials satisfy combinatorial reciprocity theorems, i.e., they have an (a priori nonob- 
vious) interpretation when evaluated at negative integers. The classical example is Stanley's theorem 
connecting the chromatic polynomial of a graph to acyclic orientations [Sta74]; the reciprocity theorem 
for flow polynomials is much younger and was found by Breuer and Sanyal [BS1 lb], starting with a ge- 
ometric setup not unlike that of our proof of Theorem 2. We hope to lift their methods into the world of 
simplicial complexes in the near future. 
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